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I. INTRODUCTION 

Despite our lack of a complete theory of quantum gravity, quantum fluctuations of the 
gravitational field can stili be tracked down in some approximate regimes. For this reason, 
semiclassical approximations play a key role in attempts to characterize quantum gravitational 
effects, and WKB methods were developed early on in canonical quantum gravity studies [T]. 
Such methods are not only useful to obtain corrections to classical general relativity [2], they 
are also a crucial tool for tackling conceptual problems such as the problem of time or the 
emergence of a classical background metric |5H8]. Unfortunately, Standard WKB techniques 
fail drastically near a turning point. 

Surprisingly, despite their occurrence in many scenarios, the generic structure of turning 
points in quantum gravity has not attracted the attention it deserves IMI]. But turning 
points are not exotic. They arise in many cosmological models, when the Universe undergoes a 
bounce na [13], or recollapses due to its spatial curvature na, but also in black hole evaporation 
scenarios such as in 1+1 dimensions mm- One way to remedy the problem of turning points 
is to use a different set of variables to build semiclassical Solutions. In Standard approaches the 
wave function takes the spatial metric and matter degrees of freedom as arguments. Instead of 
this metric representation, one can work in a momentum representation. Since metric variables 
are canonically conjugated to extrinsic curvature variables, the wave function will now take 
curvature degrees of freedom as arguments. The idea to use such a momentum representation 
has recently attracted more attention [T61 - IT9] . While most of these works followed a reduced 
phase space quantization procedure (the “identify time before quantization” of Isham’s classih- 
cation m), our aim is to start from the fully quantized Wheeler-DeWitt constraint and control 
the validity of the semiclassical momentum representation. In a previous work [19] we were 
concerned with the dynamics of matter helds (“light degrees of freedom”) near a turning point. 
Here instead, we focus on the properties of the gravitational part (“heavy degrees of freedom”). 
The objective is to derive the main properties and demonstrate the usefulness of the momentum 
representation in quantum cosmology. 

We shall compare these Solutions of the Wheeler-DeWitt equation, and their validity condi- 
tion, with Standard WKB Solutions, and usual resolutions of a turning point singularity. The 
momentum representation provides much interest compared to alternatives, which we derive 
and discuss carefully. To start with, it is the most semiclassical variable in the vicinity of a 
turning point, but it is also well dehned far from it. As we argue, there is no necessity to “go 
back to the geometric representation”. Moreover, it is well-known that an internal variable can 
play the role of time efficiently only if it behaves semiclassically enough (this has been obtained 
from a variety of approaches: internal correlations [20], effective techniques [21], and matter 
transitions [221 ES])- The semiclassical Solutions we study are parametrized by the momentum, 
and therefore, they allow us to use it as a time variable. We apply this to several examples, 
and in particular, we show how to apply our framework to the York construction of “extrinsic 
time.” York time played a key role in the demonstration of the well posedness of the Cauchy 
problem in general relativity [21125]. It is also a central ingredient in some recent modihed 
gravity theories [20], and was recently used in cosmological settings m EH]. Here we show 
that our method exactly allows us to control the approximation necessary for York time to be 
a good time variable at the quantum level. In the last part of this work, we show that, despite 
their semiclassical character, these Solutions also allow us to access nonclassical quantities, such 
as tunneling amplitudes. 

In the hrst section, we review the construction of the (semiclassical) momentum represen¬ 
tation, together with the Standard (WKB) semiclassical approximation. In the second section, 
we discuss the key points of our construction, namely its link with the metric representation 
and the generic structure near a turning point. We then present applications of our framework. 
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In the last section, we show how tunneling amplitudes can be obtained from our framework. 


II. WHEELER-DEWITT EQUATION AND MOMENTUM REPRESENTATION 

A. Minisuperspace Wheeler-DeWitt equation 

When we restrict ourselves to cosmological systems (homogeneous and isotropic), gravita- 
tional degrees of freedom reduce to a single one. The set of possible metrics in this niinisuper- 
space sector is then given by 

( 1 ) 

where dQ‘j^ is a fixed spatial metric of constant curvature K. The coupled system of matter 
and gravity is then governed by the minisuperspace action ^ 

^ ^ / (e P) 

where idmat(a;0j) is the Hamiltonian density of matter fields, and Vu fhe volume of the Uni¬ 
verse V{a) = —Ka^ + is the Wheeler-DeWitt superpotentiah In this setup, it is well 
known that if the Universe is closed and have a nonzero cosmological constant, there is a turning 
point [iniiiiiiii]. However, a turning point may arise for a large variety of reasons. It appears 
when the Universe contains classical exotic matter fields (see e.g. [12]), or quantum matter 
with negative energy density (e.g. [2H|) or also in Loop quantum cosmology, where quantum 
effects forbid the Universe from having a volume smaller than a minimum value m- Our 
analysis aims at staying rather general, and for this reason, we shall replace the superpotential 
by an “effective” potential I4ff(a) that contains quantum effects or contributions from helds at 
equilibrium one wishes to include. 

After canonical quantization, the quantum dynamics of the system defined in ([^ is given 
by a single Hamiltonian constraint: the Wheeler-DeWitt equation, 

Again, to keep the discussion general, we consider a large class of possible matter Hamiltonians, 
containing all matter degrees of freedom, but also gravitational perturbations [29|. However, 
for pedagogical purposes, it is instructive to have in mind specihc and simple examples. For 
instance, in a Universe hlled with a single scalar field in a potential Ucai, the matter Hamiltonian 
reads 

^mat = ^(-a5 + aVscal(0)). (4) 

In this Work, we shall make the (drastic) assumption that the Universe is filled with fields at 
equilibrium. This means that we assume the wave function to factorize as |T(a)) = \I'(a)|mat), 


^ To lighten the notations, we have redefined the gravitational coupling constant G = , where Gn is the 

Standard Newton constant, and Ac = 3A the Standard cosmological constant. 

^ Vu appears as a global factor in the classical action (i). and therefore, does not affect the classical trajec- 
tories. At the quantum level, as we shall see, the volume changes the regime of validity of the semiclassical 
approximation: the bigger Universe, the more classical it is. In fact, all the error bounds we provide in this 
work scale as This fact was used in a different context to advocate for the efhciency of effective equation 
of loop quantum cosmology m- 


= 0 . 


(3) 
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where Imat) is the state of matter fields in equilibrium (we shall relax this assumption in section 


consider several massive and massless fields at equilibrium, we have approximately 

(-^mat) = Pm + (5) 

a 

Pm is the energy density of nonrelativistic matter and pr that of massless (radiation) fields. This 
applies at equilibrium, when equipartition is valid. On the contrary, during infiation, under 
slow roll conditions, the main contribution of the infiaton comes from the potential part, i.e. 
= a^l4cai(0)/2 (aud 0 ~ const), giving a cosmological constant contribution [30] . 

In the sequel, we assume a general parametrization of gravitational degrees of freedom. The 
geometry is governed by a single metric variable h, such that a = a{h) in Eq. Q In terms 
of h, the Wheeler-DeWitt equation takes the very general form 


IIC). Doing so, we can replace hfmat by its mean value (hfmat) 0 [T]- For instance, when we 


[-dl + W{h)\ *(fc) = 0, 


( 6 ) 


For this, we first notice that daada = (a^^) ^dha{a^h) ^dh- Then, we redefine the wave function 
as 


'F(a) -)■ 


T(h) 


^ya{h)/a^h{h) 

By applying this change of wave function to Eq. ([^, we obtain the pseudopotential 


^ih) = (Vesia) + 2Ga{Hm.t)) + 


^^3a^aha%h 






4G^a^a% 


Wa{h) 


W^{h) 


(7) 

( 8 ) 


We see that the pseudopotential contains two contributions. The first one is the classical pseu¬ 
dopotential Wci, written in terms of the general metric variable h. The second is the “quantum 
pseudopotential” Wq. This contribution directly comes from the noncommuting character of a 
and da- While is insensitive to a specific choice of ordering for the Wheeler-DeWitt equa¬ 
tion, Wq is directly affected. In Eq. ([^, we have chosen the Laplace-Beltrami ordering. This 
possesses the interesting property that the Wheeler-DeWitt equation is independent of the way 
we parametrized the geometry in Eq. ([^. However, there exist many other legitimate choices 
for the ordering. Again, for the sake of generality, we shall simply keep Wq unspecified, and 
therefore our results will be valid for any choice, provided that one uses the appropriate Wq. 


B. The momentum-WKB approximation 


We shall now analyze the Wheeler-DeWitt equation under its very general form ([^ so as 
to build semiclassical Solutions. The Standard procedure consists of obtaining WKB modes. 
Those are given by 


xW 


f Pci(h')dh' 

y2pd(/!) 


(9) 


^ Some Standard choices would be to use a itself, but also h = a = ln(a). Choosing a different function a(h) 
means that one uses different coordinates on the metric superspace. We shall exploit this for instance in 
Sec.lIVB 
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We denote the WKB modes by y rather than d' so as to keep in mind that x constitutes an 
approximate solution. is the (classical) conjugate momentum of h, and obeys the Hamilton- 
Jacobi equation 

Pciihf + W{h) = 0. (10) 

Eq. ([^ provides a semiclassical solution of the Wheeler-DeWitt equation ([^, and has been 
extensively studied in the literature (see lai and references therein). However, Eq. (|^ is not 
an exact solution. It provides a good approximation when 


Erh(h) 


dhPc\ 

Pcl 


< 1 . 


( 11 ) 


The error function Er^ gives a bound on the local error due to the WKB approximation ^T] . 
However, when considering a scattering problem, the accumulated error might be much smaller, 
due to destructive interference effects [221 EI]- As far as scattering is concerned, the criterion 
0 is rather conservative. In order to keep the discussion general, we shall stick to this 
criterion, which guarantees that the error is at most of the order of Erh. 

Unfortunately, the WKB approximation of Eq. (|^ breaks down in many physically relevant 
situations. In particular, it faiis dramatically near a turning point [31132], where pc\ = 0, or 
equivalently, W{h) = 0. Notice that the notion of turning points is classically independent of 
our choice of metric variable h, i.e. of the function a{h). Indeed, we see from Eq. (|^ that if Wd 
vanishes at some point, it does for any function a{h), since different HAds for different a(h)’s 
are proportional. Wq on the other hand might affect the location (or presence) of a turning 
point, although its contribution is often subdominant with respect to Wd- 

When switching to the momentum representation, turning points will no longer be singular 
point of our approximate Solutions. In the following, we shall present a slightly generalized 
version of the framework developed in [19]. In the next sections, we will study in detail its 
properties and some applications. To proceed, we take the Fourier transform of the wave 
function 

^{Ph) = j (12) 


Assuming that this Fourier transform is well dehned has nontrivial consequences First, we 
assumed that we integrate over all real values of h. This creates difficulties if one wishes to 
use a itself as the metric variable (we shall return to that point in Sec. IIIB). Second, Eq. (12) 
implies that we discard the growing mode in the classically forbidden side of the turning point. 
As shown in [16], a well-defined momentum representation necessarily implies such a selection 
rule. It is also necessary to recover the background field approximation near a turning point, see 
the discussion in Sec. IIC After the Fourier transform, the wave function obeys the Wheeler- 
DeWitt equation in the momentum representation, i.e. 


[pl + w(td,j]i(p,) = 0. 


(13) 


To obtain semiclassical Solutions for this equation, we assume a semiclassical ansatz 


'^{ph) = A{ph)exp 




hci{p)dp ) . 


(14) 


^ It is also instructive to notice that Eq. (12) (combined with Eq. 0 ) is the unitary implementation of the 
canonical transformation {a,Pa) —)■ {ph-, ^), where a = ln(a), as is described in |16) in a similar context. 
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The key point is that this expression gives a good approximation when the amplitude is a slowly 


varying function of ph- Therefore, we shall solve Eq. (13) using a gradient expansion. The phase 


of Eq. (14) is the classical action, and is obtained from the Hamilton-Jacobi equation 


Ph + W{hAPh)) = 0. 


(15) 


This is the same equation as (10), except that we now solve for h as a function of ph instead 
of the converse. If the pseudopotential is monotonic, the solution hdiph) is unique. If W{h) 
vanishes, say at h = htp, there is a turning point. The region where W{h) is positive (without 
loss of generality, we assume it for h < htp) is classically forbidden. Hence, when ph runs from 
—oo to + 00 , hciiph) goes from +oo to htp and back again. 

Before using any approximation, the ansatz (14) simply defines the amplitude A{ph). When 
applied to Eq. (13), it gives 


Ph + '^{l‘ci(.Ph) + iSp„] A(Ph)^0. 


(16) 


Since the amplitude is assumed to vary slowly, we expand the operator in Eq. (16) in powers 
of idpf^. For this, as we detail in Appendix one must use a Taylor expansion of functions of 
non-commuting arguments. Here, the hrst order expansion reads 


w(h,M + tdp,) = w{h,M) + w'{KM)tdp, + '-w"{h,ip,))KM + o (ajj. (17) 


(By convention, W' = dhW and = dp^hd-) Using this, Eq. (16) becomes 


pI + W{hdiPh)) + W'{hd{ph))idp^ + 7:W"{hdiPh))KiiPh 


A{ph) = 0 . 


(18) 


When sorting these terms in gradients, we see that the first order one is nothing other than the 

tude 

(19) 


classical Hamilton-Jacobi equation (15). The second order part of Eq. (18) gives the amplitude 

A(Ph) = \W'{ha(Ph))r'''‘- 


We then deduce the semiclassical solution of the Wheeler-DeWitt equation 

Q-i f f^cliPhWh 


x{Ph) = 




( 20 ) 


Postulating (|20|) as an exact solution corresponds to the reduced phase space quantization 
approach 


, or “identify time before quantization” [3]. In other words, this is what one 
obtains if one first solves the Hamilton-Jacobi equation for and then quantizes. This is not 
what we obtain starting from the Wheeler-DeWitt equation ([^. As we show in appendix [A| 
(20) is a good approximation if we have 


Erp(p/,) = 




W'2 


< 1 , 


( 21 ) 


where W and its derivatives are evaluated at h = hci{ph)- It is remarkable that this condition 
cannot be guessed from the somewhat naive statement that the WKB approximation is valid 
when the “phase varies much more slowly than the amplitude,” which wouid lead to an analog 
of Eq. ( pT| like \W''/{hdW')\ -C 1. In addition, Eq. (21) shows that the error vanishes on the 
turning point. Of course a statement like “momentum WKB is exact on a turning point” has 
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no meaning; the error must be small for a sufficiently large interval. But in close vicinity of the 
turning point, momentum WKB is always a good approximation. 


In the above construction, the prehx 


“semi” 


of “semiclassical” should not be taken lightly. 
Indeed, semiclassical Solutions of Eq. (20) stili encode many quantum features of the gravita- 
tional held. In particular, no classical background metric exists at the level of Eq. (20). In 
fact, the validity of our construction, governed by Eq. (21), is independent of internal (i.e., 
matter helds) degrees of freedom. As we shall now see, this is not the case when considering 
the classical background limit. In order to recover a classical space-time, we shall use matter 
helds as probes. They will experience a classical background as long as energy changes are 
small. 


C. The background as perceived by matter fields 


To physically interpret the semiclassical Solutions obtained in Eq. (20), we use a matter held 
as a probe. In other words we ask, how do matter helds perceive the wave function of the 
gravitational degrees of freedom? For this we consider a massive held, i.e. I4cai(0) = Mcf)^ in 
Eq. To simplify the discussion, we assume that this held evolves adiabatically, that is, no 
particle production occurs while the Universe expands. This might be a bad approximation in 
realistic scenarios. However, it drastically simplihes the discussion, without altering the main 
conclusions. In [19], we considered particle production due to interactions and the conclusions 
we shall draw are maintained. In the adiabatic limit, matter states are adequately described 
by their decomposition in the adiabatic Fock basis [33l [3l| . To build this basis, we consider the 
instantaneous eigenvectors of the matter Hamiltonian 




( 22 ) 


The quantum number n is conserved in the adiabatic limit. Hence, the Wheeler-DeWitt equa- 
tion (j^ decouples into several second order ordinary diherential equations, one for each value 
of n. We then apply all the preceding results using several pseudopotentials 


vr„(/i) = -aW(v'(o) + 2GaMn) + W^(h). 


(23) 


For each pseudopotential, we construet the semiclassical solution Xn{Ph) using Eq. (20) and 
the corresponding Hamilton-Jacobi solution hn{ph)- We then obtain the general semiclassical 
solution of the Wheeler-DeWitt equation as 




Q-i S hn{p'^)dp'^ 

Cn -- \n), 


VWl 


(24) 


where CnS are constants. Such a solution is an arbitrary superposition of matter eigenstates 
associated with a semiclassical state for the gravitational degrees of freedom. When adopting 
the point of view of matter, the interpretation of this solution is rather ciear. The semiclassical 
wave function for gravity is analogous to i.e. it describes the “time evolution,” where ph 

plays the role of time. From this point of view, there are several crucial advantages with respect 
to Standard discussions of the WKB interpretation of quantum cosmology H EI ES]. First, 
since we obtained semiclassical Solutions in the momentum representation, they are perfectly 
adequate to describe physies near a turning point (as shown by Eq. (HU)- Second, we can 
consider not only one semiclassical solution, but a superposition of several of them. In other 
words, Eq. (24) does not describe matter in a single background. Each matter state perceives 
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its own background. Third, one is not restricted to consider one matter eigenstate, but any 
superposition of eigenstates. In the presence of interactions or non-adiabaticities, quantum 
transitions occur between these matter states. Quantum transitions were discussed in [22l [23] 
in the metric representation and in [19] in the momentum representation. The notion of a single 
background metric emerges from Eq. (24) only when the spread in matter energy is small, i.e., 


AE E. As in [T9l [22l |23|, we consider energy changes at hrst order around its mean value. 
We then have 

hn{Ph) ^ h{Ph) + dEh{ph){En “ E), (25) 

where h{ph) is the solution of the Hamilton-Jacobi equation ( jlsj ) for E = E. The hrst factor in 
the second term gives rise to the background notion of time via the Hamilton-Jacobi relation 


t(Ph) = / dEHPk)<ip'h- 


(26) 


This is exactly what is needed to change the phase of Eq. (24) from f dp'^^ into J dt'. Therefore 
the solution (24) becomes 


-i f h(p'^)dp'i^+i f E{t')dt' 


Cr7 6 


-fEnit')dt' 


n). 


(27) 


This now corresponds to a superposition of matter eigenstates living in a background metric 


characterized by ph (t) from Eq. (26), where the geometry is sourced by the mean value of matter 
energy E [6|. The time t that emerges from Eq. (26) coincides with the “WKB time” 


and our procedure is perfectly valid near a turning point. Indeed, we recognize in Eq. (27) the 


Standard from of a Born-Oppenheimer approximation, where T = x{Ph)4’{4^')Ph)i he. the wave 
function of the heavy part factorizes. A key aspect of this approach is that t is really the time 
as perceived by matter helds. It arises as the conjugate momentum of matter energy We 


underline once more that no background time (or metric) exists at the level of Eq. (24). Hence 


using Ph as a time variable is more fundamental than using the WKB time f of (26). However, 


one stili needs to consider a semiclassical approximation to have a sensible notion of evolution, 
as a phase for each matter state. This is even clearer when considering quantum transitions, 
as Eq. (24) induces a unitary evolution for these transitions. In [19] it was obtained that the 
evolution is unitary when the momentum WKB approximation is valid and the various h„(p/i)’s 
describing semiclassical trajectories have the same monotonicity in ph- 

At this point, we emphasize that we do not claim that ph is a fundamentally better variable, 
or “time”. Very similar conclusions to that drawn in Sec. II Cj were obtained using WKB 
methods in metric representation [ZUSmiEal ES]. However, the validity conditions in both 
representations are quite different, as shown by comparing Eqs. 0 and ( [2T| . In particular, 
close to a turning point, ph is the most semiclassical variable, and hence gives the best time 
variable. 

We now discuss the validity of the background held approximation, i.e., of Eq. (25). Since 


the matter Hamiltonian appears linearly in the total Hamiltonian constraint, so does the matter 
energy in the pseudopotential of Eq. ([^. Using the dehnition of the Hamilton-Jacobi solution 
hci, we compute the condition for the hrst order expansion in the energy change AE: 


AEdlK 

dE^cl 


AE\W"dEW -W^OeWYI 


< 1 . 


(28) 


® Of course, even in a fixed background, there are many choices of a time coordinate. The choice comes here 
from the identification of matter energy in the Hamiltonian contraint. By defining for instance aiJmat instead 
of Hmat, one obtains the conformal time instead of the co-moving time [II1I3S]. 
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Unlike Eq. (21), this condition depends on the matter degrees of freedom through the energy 
fluctuations AE. In a large Universe, the total inertia of gravity, governed by E is mnch larger 
than energy changes, since the former is proportional to the total mass of all the particles and 
the latter of the mass of a single particle. Therefore, there is generally a hierarchy between the 
two approximations encoded in Eq. (21) and Eq. (25). The semiclassical approximation is in 
general mnch better than the (more drastic) classical backgronnd approximation. This means 
that the solntions of Eq. (20) are not only useful to recover qnantnm held theory in curved 
space-time, bnt they also encode some qnantnm corrections. 


III. METRIC REPRESENTATION AND BOUNDARY CONDITIONS 
A. Inverse Fourier transform at the saddle point approximation 


Having obtained semiclassical solntions in the momentnm representation, it is instrnctive 
to relate them to semiclassical solntions in the metric representation. For this, we start from 
x{Ph) of Eq. (20) and compnte its Fonrier transform 


T(h) = 


g-i / ^cl (Ph) dp'f^ +ihph 


dph- 


(29) 


This integral gives a (semiclassical) solntion in the metric representation, bnt one whose validity 
is controlled by Eq. (21) rather than Eq. (11). In particnlar, it is valid across a tnrning point. 


As we shall now see, when both (21) and (11) are valid, this integral coincides with usnal WKB 
Solutions as in Eq. §. To show this, we evaluate the integral using the saddle point method |37j . 
This method precisely requires 0 to be accurate. At each value of h, the value of the saddle 
point ph = p*{h) solves the equation 


h = hci(p*(h)). 


(30) 


Hence, p*(h) is nothing other than the classical momentum function p^{h) = Pci(h), solntion of 
the Hamilton-Jacobi equation (10). Moreover, the phase evaluated at the saddle point gives the 


Legendre transform of the action J hci{p'ffjdp'f^] that is, the action in the metric representation. 
Indeed, from a change of integration variable p'j^ = Pci{h') and an integration by parts follow 
the identity 

rPc\(h) 


hc\{p'h)dp'h + hpc\{h) = / pd{h')dh'. 


The resuit of the integration (29) at the saddle point approximation then gives 


T(h) = 


Pc\[h')dh' 

UUPM’ 


(31) 


(32) 


where the irrelevant global phase has been discarded. Moreover, deriving Eq. (15) with respect 


to Ph shows that h'^]W' = —2ph- Therefore, the prefactor arising from the saddle point approx¬ 
imation is exactly what is needed to obtain the WKB amplitude in metric representation and 


Eq. (32) simply becomes T(h) = x{h). This is a particnlar case of the fact that different Dirae 
represenations that are classically related by a canonical transformat ion are semiclassically 
unitarily equivalent [161 EH]- Our method allows us to control the validity of this equivalence, 
which requires both Eq. (0 and (|^. 
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When there is a turning point, there are two Solutions of the Hamilton-Jacobi equation (30) 
and therefore one must sum over these two saddle point contributions. Moreover, the Hamilton- 
Jacobi equation (15), equivalent to (30), is invariant under pc\ —^ —Pci- Hence if Pc\{h) is one 
solution, the other is —Pci(h). The resuit for the wave function becomes 


T(h) 


g* /'* p^i{h')dh'-i^ g-i /'* p^i{h')dh'+i^ 

yi2pdMl Vl2pd(fc)l 


(33) 


One solution represents a contracting Universe, the other an expanding one. Of course, one 
can only obtain this expression far away from the turning point. Indeed, the validity condition 
for the saddle point approximation basically reduces to Eq. (0. which is valid away from the 
turning point. 


B. Behavior of the wave function near the singularity a = 0 


A Standard choice to parametrize the metric is a = ln(a) [52], which presents several advan- 
tages. First, the quantum potential of Eq. ([^ vanishes, and second, while the wave function is 
only defined for positive values of a, a runs from —cx) to -|-cx). However, one could be willing 
to Work directly with the scale factor a, as it allows several simplifications (such as the super- 
potential becoming polynomial). Classically, the scale factor is restricted to a > 0. Quantum 
mechanically, it is an open issue whether one would like to keep this constraint or not [3] . This 
question is the minisuperspace manifestation of a more general issue in quantum gravity, and is 
referred to as the spatial metric reconstruction problem [5]. When constructing the momentum 
representation in Eq. (12), we had to consider all real values of a and this question appears to 


be important for our construction. 

The main problem of considering negative values of a is an interpretational one. What 


meaning should one give to a < 0? As emphasized in Sec. IIC , we adopt here the point of view 
of matter. In particular, we are not trying to answer a question such as “what is the probability 
of a having such value?” Hence, a should not necessarily be interpreted as part of a metric at 
the fundamental level, but only in the regime where matter fields approximately propagate as 
in a classical background. Since we considered scenarios with a turning point, we always have 
a > 0 in the semiclassical limit. To further support this, we evaluate the semiclassical wave 


function obtained in the momentum representation at a = 0. From Eq. (29), it follows that 


^'a=0 = 


J a.d{p')dp' 
^/27r\W'\ 


-dp. 


(34) 


Because a = 0 is within the classically forbidden region, the saddle point p* is complex. As a 
consequence, the value of the wave function is exponentially smalh 




a=01 


j-ImU* a^i(p')dp' 

27,\W‘,\ ' 


(35) 


This means that whether we want to impose that T only has support on a > 0 or not, the 
difference between the two choices is exponentially suppressed. Of course, the above equation 
is valid only if the turning point lies “far enough” from a = 0, more precisely, if Eq. 0 is 
valid near a = 0. If this is not the case, then imposing boundary conditions such as Ta=o = 0 
will deform the wave function and make the momentum semiclassical solution of Eq. ( [20| a bad 
approximation. In such a case, the most reasonable method is presumably to use a different 
metric variable such as a. 
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C. Turning point vicinity 


In a small neighborhood of a bounce, the pseudopotential is well approximated by a linear 
function 

W{h) = -K{h-W^), (36) 

where the bounce occurs at h = With this pseudopotential, the Wheeler-DeWitt equation 
(13) becomes 

[pI - indp^ + nhtp] ^{ph) = 0. (37) 

Since it involves only a hrst order derivative in p^, the momentum-WKB approximation is 


exact. The solution of Eq. (37) is given by 

^{Ph) = 


exp 


-*(p|/(3't) + htpPh)) , 


(38) 


where the normalization constant is chosen to be the same as in Eq. (|20|. It is now instructive to 


use the procedure of Sec.j^and relate this solution to WKB Solutions in metric representation. 

tion in the metric represenl 

'h(h) = / exp {-i{pl/{3K) + htpPh) + iphh) 


From Eq. (29), we write the wave function in the metric representation as 

dph 


This is the integral representation of an Airy function 


V2 

\ Hence, 


TTK, 


(39) 


^(ft) = (-''“''“(ft - Atp)) ■ 


(40) 


On the right side {h — htp > 0), when going far from the turning point, the Airy function 
reduces to a sum of oscillatory terms: 


T(h) 




6 3 


3— 


+ 


V2«:i/2|h-V|V2 ^2«:V2|/,_/,^p|i/2' 


(41) 


We recognize the sum of two semiclassical Solutions, obtained from Eq. ([^. As in Eq. (33), 
one describes a contracting Universe and the other an expanding one. On the other side of the 
turning point, for h — htp < 0, the wave function decays exponentially. This corresponds to 
the classically forbidden region, i.e., the Universe undergoes a bounce and hence never reaches 
very small densities. Note that in the scenario of a recollapsing Universe, the discussion above 
is the same modulo the replacement W —)■ — fU, that is, the classically forbidden region lies 
at large volumes. What we learn here is that close to the turning point, the wave function 
cannot be approximated by semiclassical Solutions in the metric representation. Far from the 
turning point, when \h — htp\ > Eq. (41) becomes valid and one can interpret the solution 
as a superposition of two semiclassical Solutions, corresponds to the “region of validity of 
WKB” of i. ^) 


® It is noticeable that in m, we also found that the matter part of the wave function is described by an Airy 
function in the vicinity of the turning point. However, the two Airy functions are very different. Here, it 
describes the wave function of the geometric degrees of freedom, while in [TH], it described the probability 
amplitude of creating particles close to a turning point. In particular, both define a region of “close vicinity 
of the turning point” (see the discussion after Eq. (41) here and after (44) in [H]), but the characteristic sizes 
are different. 

^ Our solution (40) corresponds to the function v\ of reference [5], and described in their section 2.3, but it 
was not identified as an Airy function. 
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The resolutiori of the singular character of the WKB approximation near a turning point 
with an Airy function, as in Eq. (40), is a Standard method, which can be generalized using 
the Green-Liouville approach mm- We showed here how the momentum representation 
is directly related to it. However, in many cases, it is advantageous to stay in momentum 
representation, and not make use of Airy functions or the geometric representation. First, the 
momentum representation is not restricted to a linear potential, as shown in Sec. II B[ unlike 
Eq. (40). Second, the solution in the momentum representation (38) is valid everywhere, unlike 
Eq. (41). As discussed in Sec. IIC one can then use this representation to follow the evolution 
of matter states all along, i.e., use ph as a time variable. This comparison with the Airy function 
also illustrates the limitations of the proposal to extend the WKB notion of time beyond the 
condition (11). In this approach [521, one uses the phase of the wave function as a time. But 
Eq. (40) is purely real, and no such construction is possible in this case. On the contrary, in 
the momentum representation, the wave function is semiclassical and ph is a perfectly viable 
time variable. 


IV. APPLICATIONS 


A. An exactly solvable example 


We now apply the preceding results in several simple examples. This illustrates how the 
procedure is implemented, and what are the ingredients that make the semiclassical solution 
(20) valid. We consider a Universe with a positive cosmological constant and hlled with non 
relativistic quantum matter at equilibrium. The energy density of the quantum matter is 
assumed to be negative. This induces a repulsive gravitational force that will generate a bounce. 
Using the metric variable h = a = ln(a), the potential then becomes 

K? 


‘r(a) = -^(Ae®”-2Gp,e’“). 


(42) 


where > 0 so that —pg is the energy density of quantum matter. To build semiclassical 
Solutions, we hrst solve the Hamilton-Jacobi equation ([T^. We obtain 


«ci(p«) = 3 


Gp. 

A 


g , ^^Pg 
A 


1 + 


M 

vdPl 


(43) 


We see that «ci decreases from +cx) to dtp = In {2Gpq/A) /3 where the bounce occurs, and then 
increases back to +cxd. We then directly obtain the semiclassical solution of the Wheeler-DeWitt 
equation using our momentum-WKB expression ([20|). Note that with the pseudopotential (42), 


the Hamilton-Jacobi equation (15) possesses an inhnite number of Solutions. However, we easily 


see from the variations of W that there is a unique solution such that Q!ci(Pq) is real. Similarly, 
at the quantum level, Eq. ( [l^ possesses derivatives in pa at all orders. However, by using a 
semiclassical treatment, we basically neglect all other Solutions except the one corresponding 


to the unique real solution of the Hamilton-Jacobi equation (15). 


We now discuss the validity of the semiclassical approximation. For this we use the error 
function of Eq. (21) with the pseudopotential (42). It is easy to see that Eip^ph) —)■ 0 for 
both Pa —?• 0 and pa —t oo. Therefore, the semiclassical solution Eq. (20) is a very good 


approximation both near the bounce and far from it. The error thus grows from the vicinity of 
the bounce to a maximum value and then decreases to zero. A quick computation shows that 
the maximum error is, up to a numerical factor of order 1, given by 

Ai/2 


Er 


Pmax 


UiP,’ 


(44) 
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see Fig. On can also consider the fluctuations of matter energy 5p around the mean value 
Pq. Using Eq. (28), we see that the background field approximat ion is valid all along if 


5pK 

Pq 


< 1 . 


(45) 


We point out that, modulo the discussion of Sec. |IIIB 
when using a instead of a as a geometric variable. 


the same error bounds are obtained 



P 


Figure 1: Plot of the error function Erp/Erp^^^ as a function of P 


VvPq 


B. Cosmology with York time 


In canonical general relativity, one can define a “canonical ” notion of time using the trace 
of the extrinsic curvature [21125]. The obtained time parameter, known as York time, is the 
conjugate momentum to the volume element In minisuperspace, and hence we 

apply our results using h = a?, ph giving York time. Indeed, using the action Eq. (|^ written 
in terms of h = a^, one hnds 


Fu a 
3GYa’ 


(46) 


and therefore, ph = FjT. From the results of the preceding sections, this means that at the 
quantum level, York time becomes a good time variable only if Eq. (21) is satisfied. From 
Eq. ds]), we deduce the pseudopotential 


W{h) 


ii _ JL , iih , 

9^2 V h ) 4G2h2' 


(47) 


A peculiarity of cosmology with York time is that the pseudopotential becomes flat for large 
h. Hence, T will vary slowly and one expects the semiclassical approximation (20) to become 
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inaccurate. This can be confirmed using the validity condition of Eq. (21). York time might be 
a good time variable to follow the evolution of the state through the turning point, but it will 
eventually become a bad time variable. On the contrary, the volume, or h, will asymptotically 
be better and better, since WKB in the metric representation becomes exact for h —)■ oo. 
Therefore, an accurate description will in this case involve a mix between representations. To 
illustrate this, we investigate again the preceding example, i.e. a cosmological constant and 
Pm = —Pq (and neglecting ITq oc h~‘^). Once again, the Hamilton-Jacobi equation (15) is fairly 
simple to solve, and we obtain 

2Gp, 


hci(T) = 


(48) 


A-(3GT)2' 

As in the preceding example, the Universe contracts, bounces at T = 0 and expands again. 
However, we see here that —)■ oo at a Enite value of T (a behavior that was previously 
noticed in [IH])- This is simply due to the fact that the extrinsic curvature is given by a/a, 
which is constant in a De Sitter space-time. This means that the Universe reaches a state of 
infinite volume in a finite York time. However, this should not be interpreted literally. As one 
can see from the results of section IIC any matter field will experience its proper time, at least 
in the background field limit. In other words, it will oscillate an infinite number of times before 
h reaches oo. 

We now analyze the validity regime of the semiclassical approximation. Near the turning 
point {hci = 2Gpq/A), the momentum representation is accurate. When increasing T^, Erp 
grows until the approximation breaks down. On the contrary, Erh is very large close to the 
turning point, but decreases with increasing values of T^. At some value Ter, or equivalently, 
hcT = hciiTcr), the two error bounds cross (see Fig. [^. On one side of this Crossing, the 
momentum representation offers the best approximation, on the other side, the metric one is 
more accurate. 



Figure 2: Plot of the error functions Erp/Ermax and Erh/Ermax as a function of P = GT/A^/^. 
As in Fig. the functions are symmetric about P = 0, so we only plotted P > 0. 


To obtain a globally defined solution, valid for all values of T, we construet a combination 
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of both representations, so as to exploit the most accurate one in each range of T. For this, we 
build 


§(r) = ©(T^ 


X{h)e-^^^ 



+ Q{Tl 


T^MT), 


(49) 


where x is given by Eq. ( [20| ) and y by Eq. (§. Now two remarks are in order. First, by 
construction, the error made in this approximate solution is bounded by the minimum of the 
functions ( [II| ) and (21). A short computation shows that Ermax oc When Er^ax ^ 1, 


error 


Eq. (|49|) is a good approximatio n ever ywhere, i.e. close to the turning point and asymptotically. 

we conclude that, up to O(Ermax), 'h is smooth, and the 


IIC 


Second, from the results of Sec. 
precise value of does not matter. Indeed, both expressions coincide when both error terms 
are small, which is the case in a vicinity of T^. 

It is noticeable that the maximum error obtained using Eq. (49) coincides with the one of 
Sec. IV A This suggests that there exists a more general framework, allowing us to obtain 
semiclassical states in a representation independent manner, i.e. irrespectively of a{h) or mo¬ 
mentum vs. metric representation. The natural direction to obtain such a framework wouid 
be to use a covariant approach, i.e. path integral methods. Unfortunately, if semiclassical 
Solutions can be obtained via saddle point methods, to obtain precise characterizations of their 
validity is a very difhcult task. In molecular physics, where similar questions arise, this point 
seems to remain unclear (see na and the discussion in Sec. IV A of mi)- 


C. Tunneling amplitude from the momentum representation 


As we have seen in Sec. |IV A it is not necessary to go back to the geometric representation. 
Depending on the problem at hand, one might want to stay in the momentum representation 
all along, so as to use ph as a time variable to parametrize the rate of events (such as particle 
production [H]). However, in certain cases, going back to the geometric representation from 
the momentum allows one to go beyond the semiclassical approximation. This is the case when 
there are two distinet classically forbidden regions, and the Universe can tunnel from one to 
another. As we now show, the integral expression obtained from the momentum representation. 


i.e. 


Eq. (29), naturally gives the tunneling amplitude. To illustrate this, we shall investigate 


in detail the simplest potential barrier: the harmonic one. We consider the Wheeler-DeWitt 
equation (j^ with the potential 

W{h) = lU^ax - Xh\ (50) 


Using this expression in the Hamilton-Jacobi equation (15), the Solutions read 


hciiPh) = ± 


Wrr 


+ pI 


X 


(51) 


By picking a sign, one selects one semiclassical branch, living either on the left or right side of 
the turning points h = ±htp = i-^/WnaxT^- Without loss of generality, we shall choose the left 
branch. As we shall see, the Fourier transform of momentum semiclassical Solutions gives us 
the amplitude to tunnel on the other side. From Eq. (51) we directly deduce the momentum 


semiclassical solution (20) 


^{Ph) = 


exp ( V^max+P^dp) 


Ww^^ + pD) 


1/4 


(52) 
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where po is a conventional reference momentum value, chosen to be 0 for convenience. We can 
now build the inverse Fourier transform as a solution in geometric representation 


^{h) = 


exp 


+ *a-v2 


Ww^.. + pI)) 


1/4 


dph- 


(53) 


We directly see that when h < —/itp, we may use the saddle point approximation directly as in 
Sec. III A, and obtain a semiclassical solution reflected on the potential barrier (as in Eq. (33)). 
For h > htp, the wave function in Eq. (53) is a priori nonzero. In fact, its asymptotic for 
h —)■ +CX) gives the amplitude of the wave that has tunneled across the potential barrier. 
However, this asymptotic is slightly delicate to obtain, as the saddle point equation seems to 
have no solution for h > htp. In fact, if we analytically continue hc\{ph), the saddle point p* 
lies on the other side of the hranch cut of the square root. Therefore, to pick it up, one must 
hrst deform the contour across the branch cut so as to go to the second Riemann sheet, where 
the square root flips signs (see Fig. [^. Hence, we use the identity 


'0i 


VW^max+P^dp = 


rp* 


+ p^dp- / VlFmax +P^dp, 


(54) 


'02 


where the contour C is dehned in Fig. We can now use the saddle point method to evaluate 
(53) and we hnd 


'h(h > htp) ~ At- 


o} f Pd{h')dh' 

yi2pd(/i)i ’ 


(55) 


where Pci(h) = — IFmax- This describes the part of the wave function that has tunneled 

across the potential barrier. At is the tunneling amplitude given by 




(56) 


The tunneling amplitude is thus governed by the contour integral /g ^IFmax + P^dp. A short 
calculation then shows 


^IFmax+P^dp = 2IFmax / V^^dt, 


Hence, the tunneling amplitude reads 


JO 

.TT 

~ ^ ^IFmax- 


(57a) 

(57b) 


2 _ „-7rWmax/VA 


\At\ — e 


(58) 


This is the same resuit as one would have obtained by working directly in the geometric repre¬ 
sentation and using connection formulas. More generally, the fact that analytic continuations 
of semiclassical Solutions give tunneling amplitudes is a well-known feature of semiclassical 
methods [12] , but to our knowledge, it has never been obtained from semiclassical Solutions in 
momentum representation. This provides an alternative derivation. The interest is that one can 
keep working in momentum representation to investigate “local” phenomena, such as particle 
creation near the turning point (T^j, but the solution will stili keep track of the tunneling across 
the potential barrier, and the amount of tunneling can be obtained in fine. 
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Figure 3: Contours in the complex plane of ph- On the left: tlie real line is deformed into tlie 
represented contour, in order to evaluate (53) witli the saddle point method. The new contour 
goes across the branch cut to pick up the saddle point p*. On the right: contour going from 0 
to 0 on the other Riemann sheet, and used to obtain the tunnel amplitude in Eq. (56). 


V. CONCLUSION 


In this Work, we analyzed the general framework of [19], which gives semiclassical Solutions 
of the Wheeler-DeWitt equation in the momentum representation. The aim was to establish 
the key advantages of working in momentum representation. First, it provided Solutions whose 
domain of validity are orthogonal to the Standard WKB approximation, and in particular, are 
valid around turning points. This was shown by providing a bound for the error made by 
our semiclassical approximation and comparing it with the WKB one, see Eqs. ( [2T| and (11). 
When both approximations hold, we showed that they give equivalent results (see Sec. III A). 


Second, one needs not to go back in geometric representation. Indeed, when (21) is satished. 


the momentum semiclassical Solutions not only provide accurate Solutions, but the momentum 
becomes a good time variable, that can be used to investigate the effective dynamics of matter, 
something impossible with other turning point resolutions based on Airy functions. Third, 
despite its semiclassical character, a proper analytic continuation of these Solutions gives access 
to nonclassical quantity, such as tunneling amplitudes (see Sec. IVC). Fourth, we presented 


this method in a very general framework (for any choice of ordering of the Wheeler-DeWitt 


constraint, and any choice of parametrization of the metric, see Sec. IIB), allowing for numerous 
applications. In particular, in Sec. |IVB[ we showed how to apply our framework to the use 
of York time starting from the Wheeler-DeWitt equation. Unlike in the “choosing time before 
quantization” approach [3], where one hrst solves the Hamiltonian constraint for the desired 
variable and then quantizes, we control here under which conditions this ad hoc procedure is a 
good approximation to the general Wheeler-DeWitt equation. Since York time is the conjugate 
momentum to the volume element iza, our framework directly gives semiclassical Solutions 
parametrized by York time, and allows to control their validity regime. The study of our 
examples also hints toward the existence of a unihed framework, allowing us to minimize the 
error of the semiclassical approximation. We also believe our results are directly applicable to 
loop quantum cosmology, in the so-called Yrepresentation 
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Appendix A; Gradient expansions 
1. Noncommutative Taylor expansion 


In this appendix, we shall give a proof of the noncommntative expansion used in Sec. IIB 


(To simplify, we drop the indices used in the core of the text, i.e. Ph = P and hci{p) = h{p).) 
The aim is to expand the operator W{h{p) + idp) in powers of idp. Note that this operator 
is perfectly well defined. It can be built rigorously as the composition of the multiplication 
operator Mw : 'ip{z) — )• W{z)'ip{z), with the Fourier operator and the (unitary) multiplication 
operator Mexp : '^(z) — >■ exp(—if h(z)dz)i/^(z). We then define 


W(h(p) + tdp) = 


^exp‘ 


(Al) 


To obtain its Taylor expansion, we assume that the function W{z) is analytic, and possesses 
the following series representation 


W{z) = Y, 


ar,.z 


nSN 


The operator we wish to expand is therefore given by 

W(h(p) + idp) = (h(p) + idp)' 


(A2) 


(A3) 


tiGN 


Expanding this equation in powers of idp is a well defined procedure. By expanding the series 
term by term, we shall see that it is rather easy to sort out the various powers of idp. The 
delicate point however, is to understand the validity of a truncation of the series, as being a 
“good approximation” of the full operator. In the next subsection, we return to that point and 
characterize the error made by such a truncation using a next-to-leading order computation. 


Back to Eq. (A3), its n-th term gives 

{h{p) + idp)'' = hip)" + + oidl). 


n—1 


(A4) 


m=0 


The second term contains all the possible orderings of idp and h{p). To simplify it, we notice 
that 


n—1 


n—1 


{HpT~^idp + imh{p)^-^h'{p)) , 


m=0 


m=0 


h'{p)) , 

(A5a) 


(A5b) 


= nh{p)"‘ ^idp + -n{n — l)h(p)"' ^h'{p), 

where the well-known identity ~ l)/2 has been used. Plugging this resuit into 

Eq. (A3) and summing up, we directly obtain 


W{h{p) + idp) = W{h{p)) + W'{h{p))idp + -W"{h{p))h'{p) + 0{dl). 
This can also be written in a more symmetric way as 

W{h{p) +idp) = W{h{p)) + ]:{W'{h{p))idp + idpW'{h{p))) +0(9^). 


(A6) 


(A7) 


This form shows that the hrst order term is simply the only ordering that preserves the Her- 
mitian character of the operator W{h{p) + idp). 
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2. Next-to-leading order and error function 


By using the same method as before, a lengthy but straightforward calculation gives the 


second order expansion of Eq. (A7) 


0(5=) = \-m,W"(h{p))td, + i ((id,YW"(h{p)) + H="(ft(p))('i5p)=) + o{S, 


4 ^ ^ ^ 8 

From this, one can obtain an estimation of the error due to the semiclassical approximat ion 


(A8) 


equation (13) as 


that is, the validity condition (21). For this, we decompose a solution of the Wheeler-DeWitt 


^(p) = X(p)(l+ e(p))- (A9) 

X is the semiclassical solution (20), and e encodes the corrections. A next-to-leading order 
gradient expansion of the Wheeler-DeWitt equation (|T^ gives an equation on the error function. 

(AIO) 

where A is the semiclassical amplitude of Eq. ([20j) and W and its derivatives are evaluated on 
h = h{p). To ease the discussion, we assume W' > 0 and since A{p) = directly 

obtain e\p) in terms of W and its derivatives. 

ASW'W"W'"h'^ + QWWh" - 


W'A{p)ie'{p) + \idpW"idpA{p) + ^ + W^\tdpf) A{p) = 0, 


e{p) = i- 


(All) 


32W'3 

Although e is a local function of p, what is relevant to characterize the semiclassical approxima- 
tion is the accumulated error. Starting at a point p = Po, X always matches an exact solution 
of the Wheeler-DeWitt equation (with appropriate initial conditions). At a later point pi, the 
relative error is then bounded by W{p)\dp (see for instance [SH [37] for the Standard WKB 
approximat ion, and the appendix of [Hj in a context of higher order differential equations). If 
the variations of e' are smooth enough, this integral is bounded by the maximum value of e. In 


addition, if the derivatives of the potential are well sorted, all terms in Eq. (All) contribute to 


the same order. The numerical factors of the various terms above are thus irrelevant, and one 
can tune them to integrate e'{p) in a simple manner. Doing so, we obtain the estimate for the 
error function as 

W”^h' 


The criterion of Eq. (21) is then simply |e| 1. Note that this condition is very general. In fact 


it would stili apply if the constraint were of the form Q{p) + W{h) = 0 rather than p"^ + W{h). 
In our case, one can use the the Hamilton-Jacobi equation p^ + W{h{p)) = 0 to rewrite (A12) 
and eliminate h'. This gives alternative forms 


e{p) ~ —i 


2W 


112 


p 


~ — 


(A13) 


W'^ W'^ 

The last form presents the advantage to expressing the error function only in terms of W and 
its derivatives. 
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